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The Ginsparg- Wilson relation is extended to interacting field theories with general linear symme- 
tries. Our relation encodes the remnant of the original symmetry in terms of the blocked fields and 
guides the construction of invariant lattice actions. We apply this approach in the case of lattice 
supersymmetry. An additional constraint has to be satisfied because of the appearance of a deriva- 
tive operator in the symmetry transformations. The solution of this constraint leads to non-local 
SLAC-type derivatives. We investigate the corresponding kinetic operators on the lattice within an 
exact solution of supersymmetric quantum mechanics. These solutions - analogues of the overlap 
operator for supersymmetry - can be made local through a specific choice of the blocking kernel. 
We show that the symmetry relation allows for local lattice symmetry operators as well as local 
lattice actions. We argue that for interacting theories the lattice action is polynomial in the fields 
only under special circumstances, which is exemplified within an exact solution. 



I. INTRODUCTION 

Lattice simulations of the path integral are a power- 
ful tool to study quantum field theories, especially their 
non-perturbative properties. The first step in this pro- 
gram is to find a correct discretisation of the continuum 
action and its symmetries (and appropriate discrete ob- 
servables). In that respect a lot of experience has been 
collected in the case of gauge and chiral theories. 

There have also been efforts to simulate supersymmet- 
ric (SUSY) theories on the lattice for many years. For a 
quadratic, i.e. free theory, a supersymmetric lattice ac- 
tion can be constructed. However, one encounters great 
difficulties in finding more general lattice actions that 
are invariant under the (naively) discretised supersym- 
metry transformations. One of them can be traced back 
to the fact that SUSY transformations contain deriva- 
tives of fields and that the continuum SUSY actions are 
invariant up to total derivative terms. In general the 
corresponding terms on the lattice do not vanish because 
the Leibniz and chain rule of differentiation is violated by 
any lattice derivative operator [l| . The optimal choice in 
this respect is the SLAG derivative [2], which preserves 
the Leibniz rule in the first Brillouin zone (BZ) (see Q 
for recent works), but leads to non- localities Q. 

Simulations without a realisation of supersymmetry on 
the lattice [5| generically suffer from fine-tuning problems 
in the continuum limit. In higher SUSY theories it is pos- 
sible to realise a part of the supersymmetry, and one can 
hope that this partial realisation already ensures the cor- 
rect continuum limit. Many approaches to lattice super- 
symmetry rely on such constructions, as e.g. [6]. Another 
possibility is to use a (nonlinear) deformation of the con- 
tinuum transformations on the lattice [7|. One then still 
has to show that such lattice transformations resemble 
the SUSY transformations in the continuum limit even in 
the presence of quantum corrections to all orders of per- 
turbation theory, as done for a specific model in Q . This 



can ensure the correct continuum limit of these deformed 
transformations, but a non-perturbative argument would 
be desirable. A further attempt to keep the full SUSY 
by representing it on non-commutative objects or link 
objects Q has been criticised in [l3|. For more details 
on lattice supersymmetry we refer the reader to [ll| . The 
conclusion of all these studies seems to be that there is 
a fundamental obstacle against interacting local lattice 
theories with exact SUSY implemented in the naive way. 

This situation resembles very much that of chiral the- 
ories on the lattice. The Nielsen-Ninomiya theorem [IJ] 
forbids exact chiral symmetry on the lattice under a few 
very reasonable assumptions like locality. Ginsparg and 
Wilson have derived a modified symmetry relation for a 
chiral lattice theory [13|; in particular, the lattice Dirac 
operator anticommutes with 75 except for a term that 
vanishes in the continuum limit a — *■ 0. A solution to 
this relation was found by Neuberger jl4|, which is not 
ultra- local [ig] , but local if the background gauge field is 
smooth (the plaquette is close to the identity) [16]. 

The Ginsparg-Wilson (GW) relation is not a mere 
modification of the naive lattice symmetry with terms 
that formally vanish in the limit a ^ 0; its justification is 
the analysis of the Wilsonian renormalisation group [I3| , 
for the general setting see also [l^ [l^ [20, |2l| . Thus it is 
a non-perturbative construction of a deformed symmetry 
transformation. Our strategy in putting supersymmetry 
on the lattice is therefore to revisit the corresponding 
procedure (this was already suggested in [22] and some 
earlier but incomplete attempts can be found in [23|). It 
leads to a lattice theory in an effective Wilsonian sense. 
This can be viewed as integrating out quantum fluctua- 
tions up to a finite number of lattice degrees of freedom, 
the integration over which is performed numerically. It is 
well-known that in such a process of quantisation classical 
symmetries get deformed. The aim of the present work is 
not only to extend the Ginsparg-Wilson approach to su- 
persymmetry. We will also generalise the investigations 
of Ginsparg and Wilson for an arbitrary linear symmetry. 



The reduction of the degrees of freedom is done with an 
appropriate blocking transformation and the results will 
depend on the choice of a blocking kernel. 

If only a quadratic theory is considered the effective 
lattice action for a given continuum action can be calcu- 
lated. This was done in [2J| for a supersymmetric model. 
Here we are not attempting to find such explicit solu- 
tions, since this is possible only for a free theory. Rather 
we investigate the implications of a continuum symmetry 
for the effective lattice action. If these implications are 
fulfilled a symmetric theory should be approached in the 
continuum limit. 

As our main result, we derive an exact relation that a 
lattice theory has to obey in order to represent the contin- 
uum symmetry. Also interacting theories are included in 
this generalisation of the GW relation to a general linear 
symmetry. By contrast, the GW relation deals only with 
the chiral symmetry acting on the quadratic fermion part 
of the action; the gauge fields are mere spectators. As in 
the GW case the naive lattice symmetry will be modified 
by terms that are proportional to the inverse blocking 
kernel. If this relation should represent a proper lattice 
symmetry certain conditions, especially locality, must be 
satisfied. These conditions exclude e.g. the Wilson oper- 
ator as a solution for chiral symmetry. 

Applying the formalism to supersymmetry with its 
derivative transformations, an additional constraint has 
to be fulfilled in order to derive a lattice remnant of the 
continuum symmetry. The solution of this constraint 
generically leads to the non-local SLAG operator in the 
lattice transformations. We also investigate the contin- 
uum limit of the blocking kernel, which is less restrictive. 

We are able to solve the relation for the quadratic (ki- 
netic plus mass) sector of supersymmetric quantum me- 
chanics (SUSYQM). In the case of this one-dimensional 
toy model the relation already yields non-trivial differ- 
ence operators - analogues of the overlap operator for 
the case of supersymmetry. 

The locality of these derivative operators can be im- 
proved using the freedom in the blocking kernel. This 
is an important result of our strategy: The blocking ker- 
nel helps in achieving the desired properties of the lattice 
theory, 'at the expense of introducing a rhs in the lattice 
symmetry relation rendering it different from the naive 
one. 

Concerning interacting SUSY theories, the symmetry 
relation generically couples different powers of the fields 
in the action beyond second order. It is therefore 
very intricate to truncate the interaction in the power 
of fields. As an example displaying these difficulties 
we solve the case of constant fields in SUSYQM. For 
general theories we give a necessary criterion for the 
construction of polynomial interactions. 

The paper is organised as follows: First we introduce 
the blocking procedure in detail and derive the symme- 
try relation. It is particularly simple when a quadratic 
action is considered, and for chiral symmetry we recover 



the GW relation. The next section is devoted to the 
additional constraint and its solutions. Their continuum 
limit as well as the continuum limit of the blocking kernel 
is analysed carefully in Sect. IIVI 

In Sect. |V]we start to apply the approach to SUSYQM, 
giving the general solution for the quadratic case and dis- 
cussing its properties, especially locality. Afterwards in 
Sect. IVII we solve the constant field case and investigate 
the possibility of polynomial actions fulfilling the sym- 
metry relation, both in general and for SUSY. We end 
with a summary of our results and a discussion of their 
implications. 



II. WILSONIAN EFFECTIVE ACTION 

In this chapter the Wilsonian effective action com- 
puted from a general blocking transformation is intro- 
duced. The key result of the following investigation are 
the blocked symmetry transformations. These are speci- 
fied later for the case of supersymmetry. The simple form 
of the symmetry transformations in case of a quadratic 
action is derived in section IIIDI With this result the 
general concepts are elucidated at the example of chiral 
symmetry where the well-known Ginsparg- Wilson rela- 
tion is reproduced. 



A. Blocking 

The starting point is the generating functional of a 
general quantum theory 
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and fields ip{x) — 
both, bosonic and 



with classical action S'ci[93] 
{tf^{x) , ...if [x)) comprising 
fermionic degrees of freedom. 

The theory is reduced to a finite number of degrees of 
freedom by introducing averaged or blocked fields 4)f {an) 
at lattices sites n E N1XN2X . . .x Nd by a linear blocking 
procedure. 



(j>f{an) := d x f{an — x) ip{x) . 



(2) 



where a is the lattice distance, and a blocking function / 
peaked at 0. / should have the dimension inverse to the 
d-dimensional integral, such that the original and blocked 
fields have the same dimension. 

This blocking can be easily generalised to finite block- 
ing steps and fields (/?* where i comprises internal indices 
and species of fields. Then a general blocking reads 
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where a summation/integration over indices is under- 
stood. The blocking matrix f^^ is rectangular, and usu- 
ally relates only the same species, but in general it also 



mixes the internal indices. In ([3]) a; can be a discrete 
index, then ip^ is aheady understood as a blocked con- 
tinuum field ipx = {4>fo)x- In that case the blocking maps 
a lattice with a smaller lattice spacing and consequently 
more degrees of freedom onto a smaller one with a larger 
lattice spacing. 

The blocking procedure amounts to rewriting the gen- 
erating functional Z in eqn. ([T]) in terms of a path integral 
on the lattice defined by the image of /. For the sake of 
simplicity we concentrate on the generating functional at 
vanishing sources for now, the general case, important for 
the discussion of the continuum limit, will be discussed 
in section HVl We proceed with rewriting the generating 
functional Z[0] as 
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where the 0^ live on the lattice given by the blocking /, 
and 5* is the Wilsonian effective action with 
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SDet is the super-determinant, i.e. the determinant for 
bosons and its inverse for fermions^. In ([5]) we have in- 
troduced a quadratic smearing with a blocking kernel 

{c^-4,j)a{ct>-ct>f) = (0-0/)>jf™(0-0/)j„- (6) 

By inserting ^ into ^ and performing the Gaufiian 
integration over it can be straightforwardly checked 
that (HI gives Z[0] in ([T]). Note that the above smearing 
also encompasses a's with diverging or vanishing deter- 
minants. This and further details will be discussed in 
section IIVI The explicit examples of quadratic actions 
and chiral symmetry are shown in section IIIDI and sec- 
tion III El respectively. 

For illustration of the smearing procedure we briefly 
discuss a specifically simple case with ajf^ — aS^^Smn 
with a — > cxD. Then the smearing term turns into a 5- 
function in field space. 
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and the smearing is removed. Removing the smearing 
in the continuum limit, a ^ 0, is necessary in order to 
recover the original action in this limit, S — > Scu apart 
from the blocking / becoming the delta-distribution, that 
is 0/ -^ (p. Furthermore, for achieving Q, a/a'^ has 
to diverge. More general restrictions for the continuum 
limit will be addressed in detail in section [TVl 

The blocking kernel shall connect bosons and fermions 
only among themselves and in these subspaces it obeys 
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(8) 



^ Strictly speaking SDet 'a in our case means the PfafEan of a 
for fermionic fields and det~^' a for the bosonic part. 



where 



m = 



fermionic 
bosonic 



(9) 



In other words, a is antisymmetric for fermions and 
symmetric for bosons, a = zta^ , where the minus sign 
applies whenever fermionic indices are interchanged in 
the transposition of the matrix. 



B. Symmetries 

The primary concern of this construction is the block- 
ing transformation of symmetries of the classical action. 
We now investigate in what form the lattice action inher- 
its this symmetry. Our main result will be the relation 
(fT5)l that corresponds to the Ginsparg- Wilson relation, 
but is valid for a general linear symmetry. Continuum 
implementations of the related ideas close to the present 
line of arguments can be found in e.g. [20l . |21| , for reviews 
see [11 [g. 

Let the classical action S be invariant under a linear 
transformation 
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p + Sp, {Sp>)l = eMgip: 
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where M in general relates different field species (i,j), 
but may also act non-trivially on the coordinates {x,y). 
e is the small parameter of the transformation. In the 
application to SUSY, e is Grassmann-valued as M mixes 
bosons and fermions, and M also contains derivatives. 
We have introduced the notation that for a given lattice 
quantity O the O refers to the corresponding continuum 
quantity. 

In combination with the averaging function /, see eqn.s 
(P|) and ([3]) , this symmetry transformation induces a cor- 
responding transformation on the blocked field 0^ 
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Of course we want to represent this transformation solely 
on the lattice fields. Indeed, the transformation can be 
lifted to 6 as 



with a lattice transformation A/, if 



ik fkj 






nmJ mx 



pk j^kj 
J ny yx ' 



(12) 



(13) 



holds. 

This property can be viewed as a constraint ensuring 
the compatibility of the lattice symmetry transformation 
with the blocking. One might be tempted to use it to 
define M, but / has no right-inverse, since it maps onto 
fewer degrees of freedom. This constraint has been men- 
tioned without further investigation in [23| . We will anal- 
yse it in full detail in Section IIIII and argue that it has 



severe consequences, if M contains a derivative. Since the 
transformations defined by M act on lattice fields they 
can be regarded as a naive realisation of the symmetry 
transformations on the lattice. 

According to eq. (|12p this naive transformation 5 has 
the operator representation 



<5 = eM^„0^-- 



(14) 



on the space of fields 0. In the fermionic sector the left 
derivative is used, e.g. ^V'V' = ~i^- 

Now we can transform the classical symmetry into a 
relation of the effective theory on the lattice. To that 
end we apply in equation ([5|) the naive symmetry trans- 
formation — > + (5(/) defined above to the Wilsonian 
action as well as the classical symmetry transformation 
Lp ^> Lp -\- 5(p io the integration variable using that the 
invariance of the classical action, (55ci[(/3] — 0. To linear 
order in e we arrive at 



expansion of the Jacobi determinant and comprises the 
possible anomaly of the symmetry transformation 5. 

The difference {(j) ~ (j)f) in this equation can be ex- 
pressed as a 0-derivative using 
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When inserted into (|T5|) , the first term on the rhs of this 
identity contracts to the supertrace +e STr M of the lat- 
tice symmetry A/, while the rest only contains derivatives 
w.r.t. the blocked field and can be pulled outside the 
Lp integral. Then (jlSp turns into 
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The supertrace term STr M on the rhs comes from the 
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Finally, performing the derivatives leads to a nonlinear 
relation for the blocked action S containing up to second 
order derivatives in 6, 
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(STrAf - STrAf) , 



(18) 



This is the key relation for the Wilsonian or lattice ac- 
tion S[4>], the naive symmetry transformation Af and the 
blocking kernel a. While the Ihs. of this relation is just 
the naive symmetry variation of the action 5', the rhs 
constitutes some nontrivial modification of it, that has 
been derived in the blocking procedure. The behaviour 
of this term with respect to the continuum limit will be 
investigated in section llVl Note furthermore that eq. ([18]) 
represents the lattice version of the (modified) quantum 
master equation, see e.g. [l8, 19, 20, 21, 25). 

A few comments are in order here. The supertraces of 
Af and Af -terms in this relation carry the non-invariance 
of the measures dip and d(j) respectively, and hence com- 
prise possible (integrated) anomalies of the theory. More 
precisely, STrAf carries the full anomaly related to the 
measure dip. The blocking removes a part of the integra- 
tions from the path integral leaving only the field to 
be integrated. The related part of the anomaly leads to 
STrAf - STrAf. 

The above derivation also works for a finite blocking 
step where ip is already a blocked field and x, y are lattice 
coordinates. Then, however, the starting point must be 



regarded as a Wilsonian action for a fine lattice (S'ci[<^] = 
S[ip\) that already satisfies the relation (fT5|) . which e.g. 
brings in the continuum anomaly. This again leads to 
p8)l for the blocked effective action on the coarser lattice. 
If the right hand side of p^ vanishes, we are left with 
the invariance of the action under the naive symmetry 
transformations. 
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This happens for symmetric blocking matrices as fulfill- 
ing 



Afa^^ ± (Afa 
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(20) 



since only the (anti)symmctric part of Af a~^ enters the 
rhs of the relation. The minus sign appears only if 
the matrix Af connects fermions with fermions, i.e. if 
fermionic fields are transformed into themselves by the 
symmetry. The above condition just means that the 
blocking kernel is invariant under the naive symmetry 
variation. More generally, a~^ + ag 



^ + ac;^ leads to the same 



symmetry relation ([18)) for all ag . This defines a family 
of equivalent blocking kernels a~^{ag^). 

For the chiral case a = as is excluded by the vector 
symmetry as we will elaborate in Section III El For the 
case of supersymmetry, there is in general not such an ar- 
gument and, indeed, such a matrix has been used e.g. in 
[231 • However, we shall show below, that the naive sym- 
metry M in the systematic blocking approach to SUSY 
is inherently non-local and hence excluded. Instead, a 
non-symmetric blocking kernel a must be used. 

Then the relevant symmetry can be written in terms of 
a modified field-dependent symmetry operator, Afdcf(0), 
defined as 



{Mdcf)]i„^(jji^ := Af^-'^ 0^„ - (a ^)^^, 
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Inserting this definition into the symmetry relation (jlSp . 
we are led to the relation 



(Mdef):^™0^„# = (-1) 
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the right hand side being related to a total field- 
derivative. The above derivation of A/dof follows closely 
the analogous continuum arguments as used in [l8|, [l9|, 
[20 . \2]\ . A discussion of various representations of ([Hj) 
and their use can be found in llSl. 



that Qfj^^^j is local. Locality of Afdcf reads 
lim |Afdcf(a;,y)|<e-^l^-^l 



(23) 



for some c > 
shall relax 



0. In the present investigation we 
and demand 



(24) 



|a;''A/dcf(a;,y)| < oo VreN, x,yeaN, 



for explanations see Appendix [Hj Clearly opera- 
tors Afdcf with ([23|) satisfy ([24|) but ([24|) also allows 
for softer decay, e.g. polynomial times exponential 
decay. Moreover, for interacting theories the lo- 
cality conditions P5[) .(|24 p involve field-dependent 
terms as Afdcf (0) in (|2ip is field-dependent. 



(2) Afdef has to carry the original continuum symmetry 
related to the symmetry operator M. This condi- 
tion excludes e.g. the trivial solution Afdef = 0. For 
this solution it is clear that the symmetry pattern of 
the lattice action is not entailed in Afdef, and Afdef 
does not tend towards the continuum symmetry Af 
in the continuum limit. We can summarise this 
condition in the demand that Afdef is identical to 
the continuum symmetry operator Afcont up to lat- 
tice artefacts at p = 0, where the continuum limit 
is located. Hence, the condition that Afdcf carries 
the continuum symmetry can be formulated as 



lim Afdcf 



Afco„t(ll + 0(ap)). 



(25) 



C. Local lattice symmetries 

It is important to emphasise that ([22[) in general does 
not comprise a symmetry, as the above construction ap- 
plies to any blocking kernel a. Thus, in general ([22[) 
only disguises an explicit symmetry breaking induced by 
the blocking. We shall exemplify this statement in sec- 
tion III El at the standard Wilson Dirac operator that ex- 
plicitly breaks chiral symmetry, but still satisfies ([22[) . 

The question arises what are the additional conditions 
on Afdef that make it to a deformed symmetry. For lo- 
cal continuum symmetries it is important that the corre- 
sponding lattice version of the symmetry carries this lo- 
cality. More generally, for a given symmetry the blocking 
should only induce a local symmetry breaking or defor- 
mation generated by / and a. Consequently we are led 
to two conditions: 

(1) a mandatory condition for a deformed lattice sym- 
metry is the locality of Afdcf. This guarantees a 
well-defined continuum limit in which the lattice 
artefacts related to the deformation tend to zero in 
a controlled way as they are local. Hence, in order 
to have a deformed symmetry the family of block- 
ings a~^{ag^) must contain at least one blocking 

"^ local ^h^t leads to a local symmetry operator Afdcf. 
We emphasise that this does not necessarily imply 



Note that Accent = M only for a ^ = in the con- 
tinuum, see the discussion in Section [IV Bl 



The above conditions pl|) - ([25[l should be seen as a def- 
inition of a deformed symmetry, and put constraints on 
the blocking kernel a. In order to formally obtain a sym- 
metric continuum limit, one might use actions without 
such a symmetry, but the above considerations guaran- 
tee the existence of a local lattice symmetry for every 
finite lattice spacing that converges locally towards the 
continuum symmetry. The latter property is very impor- 
tant for a successful numerical implementation. 

In the case of lattice supersymmetry the question is 
whether such a deformed symmetry operator A^dcf ac- 
cording to this definition can be constructed. 



D. Quadratic action 

The general results above simplify greatly for quadratic 
actions 



S 






(26) 



with the kernel K comprising kinetic and mass terms. At 
first sight this case seems trivial as it describes a free field 
theory. Nonetheless, it already includes the non-trivial 



case of Ginsparg- Wilson fermions [ij| with background 
gauge fields, see the next subsection. Moreover, locality 
of a symmetry operator Mdd of an interacting theory re- 
lates directly to the locality of its kinetic non-interacting 
part. 

With this action the general symmetry relation ([T8|) 
simplifies to 



(j)M^K<j) 



tr {Ma-'^) K'^ + (STrAf - STrAf) 



(27) 



In many cases the second line vanishes. In case of field 
independent transformation matrices M and kinetic op- 
erators K it anyway is just an irrelevant constant. How- 
ever, in the case of anomalous symmetries it contributes 
to the anomaly. If one considers non-quadratic actions 
the corresponding term in general becomes 0-dependent. 
The first line in (|27|) has to be valid for general fields 
6 and hence we conclude that 



m'^k±{m'^kY' 



(28) 



Again the minus signs appear on the left and right hand 
side only if fermions are transformed into fermions by the 
naive symmetry M. 

The interesting information in the symmetry relation 
is that of the propagation of symmetry breaking on the 
lattice. This propagation can be seen from 
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Af^ ± MR-^ = (a"^)'^Af^ ± Ma-^ 



(29) 



This equation high-lightens how the breaking of the sym- 
metry by the blocking matrix a and the breaking by the 
kernel K must compensate each other. It also enables us 
to read-off the general solution K, 



R- 



(30) 



Here, a^ is a general symmetry-preserving term fulfill- 
ing (pn|) . We emphasise that (|30p can also be used for de- 
termining a family a{K) for a given K. We conclude that 
pairs {K~^,a~^) are unique up to symmetry-preserving 
terms a^ . The symmetry relation can be also rewritten 
by introducing the deformed symmetry matrix A/dcf as 
defined in (I2ip . Here we find a 0- independent A/dcf with 



Afdcf := M (11 - a-^K) = -Ma, 



'K. 



(31) 



Afdef may, however, now depend on background fields via 
K and a^^, e.g. link variables if K is the Dirac operator. 
Note that ([31]) defines a family of symmetry matrices, as 
ag is a general symmetric matrix satisfying (j20p . For 

reads 



the modified symmetry the relation 

MliK±{Ml,Kf-- 







(32) 



As already mentioned in the previous section, (|32|) in 
general does not comprise a symmetry, as the above con- 
struction applies to any kinetic operator. Thus, in gen- 
eral p2p only disguises an explicit symmetry breaking 



induces by the blocking kernel. We also clearly see the 
necessity of the second condition (P5|) : for a~^ = K~^ 
we have a^^ = and hence A/dcf = 0. Then the modified 
symmetry relation p2p carries no information about the 
symmetry at hand. 

In turn, only A/dcf 's in (EB) with (|32p and the locality 
and continuum limit properties (|24p and (|25p respectively 
define deformed lattice symmetries. 



E. Chiral symmetry 

We shall first discuss the above construction and condi- 
tions at the example of the chiral symmetry. Consider an 
action out of the field multiplet of two fermionic fields: 
(j) = (■)/), i/j^). The related kinetic operator is given in 
terms of the Dirac operator 
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(33) 



and the action (pS)) reads S = a'^ 'ipVip . As explained 
above, in our units the quantities relevant for the con- 
tinuum must get additional factors of a'* to account for 
the integral. The continuum action is invariant under 
symmetry transformations generated by 
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(34) 



with 75 = 75 • Since the transformation acts only alge- 
braically on spinor indices it is easy to fulfil the constraint 
P^ . The naive transformation is just the same as the 
continuum transformation. A general blocking matrix a 
carries the fermionic anti-symmetry and reads: 



a 

nd 



-a{ 
ai 



(35) 



with a general ai. Note that in order to get a real action 
both T) and ai must be hcrmitian. Inserting the kinetic 
operator (j33p . the chiral transformation matrix (|34p . and 
the general blocking ([551) into (P7)) we are led to 



{/?,75}=I?{75,ari}/?, 



(36) 



which comes from the field dependent part of (j27p . It can 
be rewritten in terms of a deformed symmetry, cf. [22 |. 
which is according to the general definition of A/dcf in 
((3T|) given as 



A/dcf = 



75,dcf 

I..- . AT 



(37) 



(75,dot) 

with 

75,dcf = 75(l-ar^^), 75,dcf = (1 -/?aj'^)75.(38) 

In terms of the deformed 75 's the symmetry relation reads 

75,dcf/? + /?75,def = . (39) 



For hermitian aj~ and V we arrive at 75, dot = 75 def ■ 

In case of a theory with vector-symmetry the blocking 
should respect it. Hence the simplest ai is a fermionic 
mass term with mass 1/a, ai = 1/all, where 11 is diagonal 
with respect to the lattice sites and the identity in Dirac 
space. The result, 



{V , 75} = 2al?75l? , 



(40) 



is the Ginsparg-Wilson relation [131 . 

The deformed symmetry operator A/dcf from (|37p . ([55]) 
is local due to the ultra- locality of a^^ and the locality 
of v. It should, however, be noted that 75, def is not nor- 
malised, 7|def 7^ ^i a-nd even vanishes at the doublers. 
We conclude that 75, dot does not define a chiral projec- 
tion. Indeed no such normalised 75 def can be constructed 
for a single Weyl fermion, see [2g|, as a consequence of 
the Nielsen-Ninomiya no-go theorem. In the given exam- 
ple the normalisation of 75, def fails at the doublers, it is 
neither smooth nor local. 

The part of equation (P7)) . that is independent of 0, 
carries the integrated chiral anomaly. 



Tr75l? -I- (Triattico75 - Trcont75) = . 



(41) 



This constrains the continuum regularisation in terms of 
the lattice blocking and vice versa. The blocking matrix, 
((35)) . was chosen to have vector symmetry. It is thus pre- 
serving the vector symmetry of the kinetic term ([55)1 . It 
follows that Tr75a~il? = 2Tr75l? = (n+ - n_)iattico, 



where n+jg 



"latti( 



are the numbers of fermionic 



zero modes with positive and negative chirality respec- 
tively related to the (background) gauge field on the 
lattice. If we choose a blocking compatible with axial 
symmetry, the related term would vanish. Then, how- 
ever, vector symmetry is broken, and we would loose 
(background) gauge symmetry. This analysis is refiected 
in the well-known fact that Trcont75 is regularisation- 
dependent. Choosing a vector symmetric regularisation 
of the trace, e.g. 



Trcont75 := lim Trcont75e 

c— tO 



f-'Di 



(42) 



we are led to Trcont75 — n^ — n_, where n+, n_ are the 
numbers of fermionic zero modes with positive and neg- 
ative chirality respectively related to the (background) 
gauge field in the continuum. In turn, an axially sym- 
metric regularisation leads to Trcont75 = 0. The cor- 
responding trace on the lattice vanishes, Triatticc75 = 0. 
Note that the sum of zero modes on the lattice here comes 
from the functional relation involving S. In summary we 
conclude from (j4ip that full chiral symmetry in the pres- 
ence of a background gauge field is maintained iff the 
lattice gauge field permits the same difference of positive 
and negative chirality zero modes. 

As an example for an explicit breaking of chiral sym- 
metry we consider Wilson fermions with Dirac operator 



aVw = ilti sin(ap^) 4- r ^(1 - cos(ap^)) . 



(43) 



In this case chiral symmetry is explicitly broken due to 
the momentum-dependent Wilson mass. We start with 
the relation ([50)1 for general Dirac operators T). The cor- 
responding blocking kernel (cf. psp ) is given by 

-1 



-i^p-i 



'1,S ' 



(44) 



The singularity of I? ^ at the centre of the Brillouin zone 
has to be removed from a~^ g in order to guarantee the 
continuum limit of Aidcf ^ Afcont, (|25p . This is achieved 
with 

<^i!s ^ ^M ^ tr 7^P-1 + ^a-^g , (45) 

with tr7^7y = — dtJ^^, and d is the space-time dimension. 
We conclude that aj~^ is given by 



= ll-trP^i -Ac 



i,s ■ 



(46) 



with a scalar first term proportional to 11, and a sym- 
metric contribution ^ol^ g proportional to '^fj,f^{p). Note 
that the first term cannot be changed by Aa|j~g and hence 
carries the unique information about the symmetry- 
breaking part of the kinetic operator K. This part is 
the same for all members of the family a{K) of blockings 
corresponding to a given K. Restrictions on this part 
will hence constrain the class of possible lattice actions. 
For Ginsparg-Wilson fermions we have aj~ — all and 
Aoj^g = 0, that is a^^ has no symmetric part. For 

Wilson fermions the choice Aa|j"^ = leads to a non- 
local aj~ : some higher derivative of tiV^ {p) is not 
bounded at the origin and this contradicts locality, see 
Appendix [21 Furthermore this non-locality cannot be 
changed by the symmetric term Aqj^^ except in one di- 
mension. We conclude that the blocking a~^ related to 
the Wilson-Dirac operator is inherently non-local. Still a 
priori this does not entail that the corresponding Afdcf is 
non-local. However, the product of this a^^ (including 
Aa^^g) with the Wilson-Dirac operator is also inherently 
non-local and enters Mdof. Actually, the non-locality of 
Afdcf is most easily seen from the non-locality of the left 
hand side of ^^ in a Taylor expansion about p — 0. 
We conclude that there is no deformed chiral symmetry 
operator Afdcf for Wilson fermions. 

We summarise that the lattice blocking induces the 
continuum regularisation. In turn, if we have chosen a 
specific continuum regularisation, this restricts the lattice 
blocking compatible in the continuum limit. We conclude 
this analysis with the remark, that an analysis of chiral 
transformations ip ^' (l±75)/2 V' completely fixes the re- 
lations, as the related integrated anomaly is independent 
of the regularisation. This is at the heart of the lattice 
observations made in [23, Eal . 



F. Explicit solution for a quadratic action 

For a quadratic action it is also possible to solve ([Sj) 
for the effective action S explicitly. Assuming S'cif'ys] = 



8 



l,oi 



V^^xyVy^ tl^s lattice action S'[(/)] = ^(j)K(j) can be ob- 
tained via performing the Gaussian integration. It leads 
to 



K^a-afif'af + K)-'f'a 



(47) 



After some manipulations that can be found in ap- 
pendix [A] the resulting fixed point operator reads in mo- 
mentum space 



K{pk 



K{Pk+l'^) 



+ a \pk) 



(48) 
Note that such a solution of the Ginsparg- Wilson relation 
was already mentioned in [13]. It is often called perfect 
lattice action. 

In most cases f{x) is considered to be the averaging 
over one lattice spacing, e.g. in one dimension 



/(^) 



1/a if \x\ < a/2 
otherwise 

2 sin(pfca/2) 



(49) 



which means f{pk) = _^ 2l£i£ii£/£!. Such an averaging 
was applied in \2^ to construct a free supersymmetric 
(perfect) lattice theory. However, since the constraint 
()13|) was not considered there the symmetry properties 



of the resulting effective action can not be expressed in 
terms of a lattice symmetry involving only lattice fields: 
equation (jl3|) demands for the derivative operator ap- 
pearing in the supersymmetry transformations 

/ Vnm '/'(aw) = —(93(071 + a/2) — (^(art — a/2)) (50) 
•^ — ' a 

m 

and this can not be fulfilled for any V„„i since the trans- 
formation involves the continuum fields. 

To interpret the rhs of equation (j50p a new field was 
introduced in [24], which is defined to be -ip{an + a/2) 
at the lattice point an. Then the lattice fields are trans- 
formed into such fields under the supersymmetry trans- 
formations. They are, however, rather a continuum than 
a blocked lattice quantity. The correct SUSY continuum 
limit is therefore ensured in this approach because the 
lattice action is a direct solution of the blocking. But 
this property can not be expressed in terms of a lattice 
symmetry, that contains only lattice fields. A well de- 
fined lattice symmetry is, however, desirable as a guid- 
ing principle for the construction of a more general lattice 
action. 



III. ADDITIONAL CONSTRAINT 



on the symmetries M, M and the averaging function /. 
It is trivially fulfilled if the symmetry transformation 
merely acts on the multiplet indices, e.g. with 75 in the 
chiral case. 

However, whenever the symmetry transformation M 
contains a derivative - as in the case of supersymmetry 
- the constraint becomes nontrivial. The problem can 
be considered in each space-time direction separately. It 
states, that the derivative d (in M) is 'pulled through' 
the averaging function / to become a lattice derivative 
operator V (in M) that acts among the averaged fields: 



dx f {am — x) tf{x) = I dx f{an~ x) dxf{x) (52) 



for all continuum fields f(x) (neglecting internal indices 
«,j) and for all lattice points n. This constraint will 
restrict the possible lattice derivatives V to be used in 
the lattice symmetry transformations M as we show now. 
In order to satisfy hermiticity and translational invari- 
ance, V should be an antisymmetric circulant matrix 



1 



(JV-l)/2 

2^ ^ 

/=-(Ar-l)/2 



Cl ^n-r 



(53) 



with real coefficients q fulfilling c-i = —ci. For simplic- 
ity we have specialised to an odd number N of lattice 
points. The Kronecker symbol S on the rhs is periodic 
with periodicity N. 

We use a partial integration in ([5^ and a Fourier trans- 
form with a discrete momentum pq = 2'Kq/L with 5 G Z 



't'l 



^ X 

/ X . 

X • 



M ~ d^ 




i\/~ V 



A. Discussion 

In the derivation of the relation of the effective action 
there has emerged a novel constraint, eq. ITT 



Mf = fM, 



(51) 



FIG. 1: A sketch of the blocking procedure: The averaging 
function / maps from continuum fields </3 to averaged fields 
(j)f (that are connected to the lattice fields <j) via a). The ad- 
ditional constraint eq. ([13} conies about because the diagram 
of / with the continuum symmetry M ^ dx and the lattice 
symmetry M ~ V has to commute. 



and a lattice volume L 
FB]) to arrive at 



Na (for details see appendix 



f{Pq)[V(p,)-ip,] =0 VgeZ. 



(54) 



Hence, the constraint states that the averaging function 
can have non- vanishing Fourier components f{pq) for 
each wave number pq for which the difference operator 
has the 'ideal' continuum dispersion relation V(pq) = ipq. 
The latter condition means that the naive translation 



(5,e^P«")| 



Ev. 



^tpqam 



(55) 



holds for all lattice points n for this wave number. 

At this point let us stress that because of the period- 
icity V(pq+Ar) — y{pq) the square bracket in eq. (f54|) 
can vanish only once for every q mod TV, for all other q 
f{pq) has to vanish. With regard to the correct contin- 
uum limit of V we demand this to happen - if at all - in 
the first Brillouin zone, thus 



f{pq)=Q forb,|>^(l-l). 



(56) 



This limits the spatial resolution of / to | Axj ~ a which, 
however, is a natural scale in the blocking approach to a 
lattice (see (HHl) for comparison). 



B. Solutions 

Inside the first Brillouin zone the constraint (|54p intro- 
duces a kind of uncertainty relation between the averag- 
ing function / and the lattice difference operator V. For 
instance, if one demands an ultra-local operator ranging 
over only one neighbouring point, V will be proportional 
to the symmetric difference. 



ciV" 



ysymm 



1 

2a 



[Sn- 



+ l,m 



ci - sin (apq) 
a 



The dispersion relation is in this case 

V(P,) 

and the bracket in ([M]) vanishes for p 

2n/N 
sin(27r/iV) ' 



^n— 1,771/ ■ 

(57) 
(58) 



1 iff 



Cl 



(59) 



The proportionality factor ci indeed approaches 1 in 
the continuum limit and hence V approaches the con- 
tinuum derivative as V^^™™ does. As a consequence, 
/ has only the lowest (and zeroth) Fourier-components 
{pk = {-2n/L,0,27:/L}y. 



fix) = fo + fi cos{2ttx/L) . 



(60) 



Hence for this ultra-local difference operator the averag- 
ing function / is very broad as it probes the whole space 

x& [-L/2,L/2]. 



The solution allowing for the next [p = 2) Fourier com- 
ponent in / demands the difference operator to spread 
over at least nearest and next-to-nearest neighbours with 
coefficients 



^(ci(5„+i. 



On-l,m) + C2{Sn+2,rn — On-2,m 

(61) 



• IL/ft (-. 

za 
with the solutions 

27r 2 sin(47r/iV) - sin(87r/iV) 



Cl 



C2 



N sin2(47r/iV) - sin(27r/A^) sin(87r/iV) 
27r 2 sin(27r/A^) - sin(47r/A^) 

^ sm^iin/N) ~ sin(27r/iV) sin(87r/iV) 



(62) 
(63) 



which correctly approach 4/3 and —1/6 in the limit N — > 

CX). 

One can proceed in this way. The more Fourier co- 
efficients are included in the averaging function, the less 
localised gets the derivative operator. In general, V needs 
to spread to the nth neighbours to enable n non- vanishing 
Fourier coefficients in f{p) ^ . 

At the extreme, in order to make / as narrow as it 
can get, all Fourier components f{pq) (in the first Bril- 
louin zone) are needed. The constraint (|54p then leads 
to the non-local SLAC-operator by definition \2\. The 
coefficients in this case are 



Cl 



SLAC 



i-iy 



2tt/N 
sm{-Kl/N) 



(64) 



All of these operators can in principle be used inside M 
to generate the lattice supersymmetry transformations. 
But there are additional restrictions on / further reduc- 
ing these possibihties. 



IV. CONTINUUM LIMIT 

The results of the last section IIIII necessitate a care- 
ful investigation of the continuum limit. We have argued 
in Section III Al that the averaging function / needs to 
approach the delta-distribution in the continuum limit. 
Hence, approaching this limit, more and more Fourier 
components f{pq) are needed. As a consequence of the 
additional constraint ([51]) , the lattice derivative operator 
V agrees in Fourier space with the SLAC derivative for 
the increasing number of modes Pq with non-vanishing 
f{pq)- In other words, the difference operator becomes 
more and more extended over neighbouring lattice sites, 
while / gets narrower. The more neighbours are included 
in the lattice derivative V, the more demanding numeri- 
cal simulations will become. One is therefore tempted to 
use the most localised solution from this constraint. 



^ In this way derivatives interpolating between the symmetric 
derivative and the SLAC derivative are constructed. 
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Moreover, if combined with an appropriate blocking 
kernel a, the numerical effort could be reduced further. 
Thus it is advantageous to determine the general setting 
giving access to the full set of allowed blockings /'s and, 
in particular, a's. 

For that purpose we reconsider the full generating 
functional in the presence of external sources. The 
physics of the blocked fields 0/ with blocking /, ^, 
is carried by general correlation functions of this field, 
{(j)f{ani) ■ ■ ■ (f>f{anj.)). These correlation functions are 
generated by 



Zf[J] 



1 



dip e 



-SdM+j^/M 



(65) 



The correlations functions of (j)f naturally live on a lattice 
defined by n e iVi x iV2 x . . . x Nd resulting from the 
blocking /. 



FIG. 2: A sketch of the mismatch in the blocking procedure, 
if / - in order to generate a local lattice derivative in the 
constraint (|54[) - has a limited number of Fourier components: 
the averaged fields (f)f have fewer degrees of freedom then the 
number of lattice points A'' used for (j>- In other words, the 
(j>f transfer information from the continuum only to a coarser 
lattice (see text). 



A. Continuum limit of the blocking / 

In the limit (/)/ [ip] -^ if the correlation functions 
{(j)f{ani) ■ ■ ■ (j)f{anr)) tend towards the continuum corre- 
lations functions. Accordingly, the most restrictive con- 
straint coming from the comparison of the lattice ob- 
servables with their continuum counterparts is that / 
must approach the delta function in the continuum limit. 
Then the lattice theory resembles the continuum up to 
minor modifications. As stated in the previous section 
this means that equation ([&i)) must be fulfilled for an 
increasing number of modes. 

Let us work on a lattice with N points and let V be 
more localised than the SLAC derivative, i.e. have an 
'ideal' dispersion relation V(pq) = ipq up to some mo- 
mentum Pmax < 7i'(l — l/N)/a. The corresponding / has 
non- vanishing Fourier components up to this momentum. 
Doublers will appear in the spectrum of such operators; 
they can be removed as shown in App. El within our solu- 
tion for supersymmetric quantum mechanics. Analysing 
the consequences for /, however, we will argue against 
these solutions in the following. 

The momentum of the lattice fields 4>{pk) is restricted 
by this momentum cutoff of / (for an explicit formula 
see (IB.9I) V That means one introduces an additional mo- 
mentum cutoff smaller than the usual lattice cutoff. In 
other words the number of degrees of freedom (Fourier 
modes) induced from the continuum via / is smaller than 
the actual number N of lattice degrees of freedom. This 
contradicts the blocking philosophy where all the lattice 
degrees of freedom should come from blocked continuum 
degrees of freedom. To be more precise, in the defining 
equation (O, the averaged fields 0/ span a vector space 
smaller than the one of the lattice fields cj). Therefore 
some lattice fields have no counterpart (j)f. Rather, 
their contribution to the lattice action S[(l)\ is a simple 
quadratic one with kernel a. This mismatch is depicted 
in Fig. H 



Another way of stating the problem is that the blocking 
/ gives rise to a resolution (an 'effective lattice spacing') 
of 7r(l — l/-/V)/pmax > a- On this coarser lattice the 
derivative is actually again SLAC. It is very unlikely that 
it yields any improvement to work on the finer lattice 
with lattice spacing a, where not all of the degrees of 
freedom are induced by a blocking from continuum fields. 

We conclude that the lattice derivative V entering the 
lattice symmetry relation as M has to be the SLAC oper- 
ator or some degrees of freedom on the lattice will have 
no continuum counterparts. In any case relation (j54p 
must hold for an increasing number of lattice modes in 
the continuum limit if / should approach the delta dis- 
tribution. 



B. Continuum limit of the generating functional 

It is left to discuss the consequences of a the general 
choice for a. This is best done in terms of the generating 
functional Zf [J] defined in ((65l) . The path integral in ((65l) 
can be conveniently rewritten in terms of a path integral 
over lattice fields (p. In section Hi Al wc have done this al- 
ready for vanishing external currents J, and a quadratic 
blocking kernel ^{(j) — 4>f)a{(f) — <j)f), see (O, with sym- 
metry properties ([5]) and ^. This procedure is readily 
extended to the general case with non- vanishing currents 
by rewriting the source term exp J J(j)f via the quadratic 
blocking kernel. 



.J't'i _ 



g-jjQ ^/ 



/' 



-i(</.-0/)Q(</.-0/) + J0 



(66) 



Inserting ([55]) into ([55]) the generating functional Zf can 
be rewritten as a lattice generating functional 



Zf[J] 



1 



N[J] 



-s[<p]+J4> 



(67) 
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with Wilsonian action S as defined in ([5]), 



The normalisation Af[J] carries a trivial quadratic depen- 
dence on the current J and reads 



J\f[J] =JVe2 



Ja-^J 



(69) 



It reduces to J\f for vanishing current. We emphasise 
again that Zf[J] in (p7|) has no dependence on a, and 
reduces to (|4]) for vanishing current J = 0. Note also 
that the generating functional in (|67p with M[0] is the 
standard lattice generating functional. As A/'[ J] is a triv- 
ial GauBian, lattice simulations for correlation functions 
straightforwardly relate to those from (|67p . 

The above construction allows us to evaluate general 
choices of a. The blocking function / will be discussed in 
the next subsection, here we will assume / to approach 
the delta-distribution such that </>/ -^ (p. We have al- 
ready argued in section IIIAI that a diverging a ensures 
that the lattice field 4> agrees with the blocked field (pj [(p] . 
A simple a is a diagonal one, a^^ = cS^^dmn with c -^ oo, 
see also section Hi Al Then we are lead to 



SDet^^^ae- 



-4,/)a(4,-4,f) ^ j(0_0^.) 



(70) 



see (O, and the integral in ([5S)l is trivially done. More- 
over, the normalisation looses its J-dependence, Af[J] -^ 
A/", as a~^ vanishes for c -^ cx3. Finally, the Wilsonian 
action is given by 



-S[0] 



dipS{(f>-(f>f[ip])e 



-SdM 



(71) 



This can be viewed as the canonical form of the Wilsonian 
action. Note also that in this case the symmetry relation 
p8)l simplifies to the standard one, as the rhs vanishes 
with a^^ — 0. 

In consequence the limit a~^ ^ is the natural choice 
for the continuum limit. In order to see how a general a 
scales with the lattice spacing a, we rewrite the blocking 
term as 



(j>f)a{(f>-(j>f) 



(72) 



nm 



{4>-<t>f)l 



Firstly we remark that a** ^„ -^ J dfix in the continuum 
limit a — > 0. Secondly, for symmetric smearings a we 
can always diagonalise the matrix a. Thirdly, a factor a'^ 
appears because Snmlo!^ -^ ^(a^ ^ v)- We conclude that 
all eigenvalues a„ have to satisfy 



a.n/a'^ 



(73) 



for guaranteeing that the Wilsonian action tends towards 
the classical action in the continuum limit. 



In other words, the inverse a~^ has to vanish for 
S —f Sci- For practical purposes it might be advanta- 
geous to work with some vanishing eigenvalues of a^^ 
already at finite lattice spacing. Then the related eigen- 
values ck„ diverge already for finite lattice spacing, and 
the right hand side of (j70|) will be proportional to delta- 
functions for the related eigenfunctions ip„, leading to 
('0n,0) = {ipn,4'f['P])- We also note that the symmetry 
relation p8|) contains only a~^, which shows no divergen- 
cies for a„ — > CXI but zeros. Indeed we show in appendix 
FdI that (fTBjl can be derived without using a explicitly, 
and hence noninvertible a~^ are not problematic. On 
the contrary, vanishing eigenvalues of a^^ mean that the 
factor N[J] in ([69]) is actually a J- independent constant 
in that subspace of fields. 

Likewise, the eigenvalues of a^^ can have any sign 
whereas ()68|) would require positive eigenvalues of a. 

What happens if some eigenvalues of a^^ do not vanish 
in the continuum limit? As the generating functionals 
Zf[J], and in particular Z[J], do not depend on a, it is 
not mandatory that the Wilsonian action S approaches 
the classical action Sd- Assuming 4>f[(p] ^ (p, a, finite 
a in (|68p amounts to equivalence classes of actions with 
measures 



exp— 5[- 



(74) 



related by Gaussian integrations in the continuunr. On 
the level of classical actions this is nothing but the intro- 
duction of auxiliary fields </> via the equations of motion. 

The case of vanishing eigenvalues a„ (diverging a^^), 
however, is different. From the path integral in ([55)1 one 
reads off that the corresponding subspace of fields p is 
simply integrated out and the Wilsonian action S does 
not depend on the corresponding eigenfunctions tpn ■ The 
corresponding singularities of a~^ in the normalisation 
(155)1 can be avoided by considering currents J only in 
the orthogonal subspace, that is {ipn,J) = 0. Then the 
above derivations are unaltered, but with this procedure 
we have removed the (yj-modes in the singular subspace 
from our theory. A simple example for vanishing a is 
given by the Wilson mass term in a fermionic theory. 
Assume that we start with the naive lattice Dirac action 
with doublers. The blocking kernel a can be chosen such 
that it vanishes at the doublers. In turn a~^ provides 
a diverging mass for the doublers. More details on this 
and the general case with vanishing a is deferred to Ap- 
pendix FeI 

We summarise our findings as follows: in order to re- 
cover the original action Sd = 5 in the continuum limit 
f ^ 6, the blocking kernel a has to lead to a delta- 
function in field space, see ([70]) . It might, however, be 
advantageous to rely on a non-trivial classical Wilsonian 
action in the continuum limit, generated by other choices 
of a, in order to optimise the locality and hence to min- 
imise the numerical effort. 



12 



V. FREE SUPERSYMMETRY ON THE 
LATTICE 

We will now apply the blocking formalism to Super- 
symmetric quantum mechanics (SUSYQM) which is a 
supersymmetric theory in one dimension, i.e. all 'fields' 
depend only on a time xi = t. It serves as a toy model 
for supersymmetric theories. After fixing the notation, 
only a quadratic theory is considered in this section; in- 
teracting SUSY theories will be considered in the next 
section. 



A. Brief review of SUSYQM in the continuum 

The field content of SUSYQM is the multiplet 

^, = {x(i),F(i),V;(i),V^(i)} (75) 

where x ^-nd F are real bosons, '0 is a complex fermion 
(Grassmannian) and "0 its complex conjugate. (The 
length dimensions of these fields are \/L, 1/\/Z and L'', 
respectively.) 

The Euclidean action in the continuum has the follow- 
ing form 



5'ci[(/?] 



dt 



(76) 
where the first three terms are kinetic ones {F is an auxil- 
iary nondynamical field) and the last two terms represent 
a potential term for x ^nd Yukawa interactions. 

This action is invariant under the supersymmetry 
transformations 



5x = — eV' + eV' 
50 — —edtX ~ ^P 



5F ~ —edtip — edt^ 
IF. 



(77) 



H = <^dtx 
up to the following surface term 

,55ci ^jdt dt (e0(9tx + F) + ei^W{x) + e^W{x)) ■ (78) 
According to our general notation we write 



where 



M 









1 








-dt 


'dt 


-1 















Sip = (eM 4- eM)ip 



M 



(79) 









-1 








"dt 











dt 


-1 






(80) 



B. Transformations on the lattice 



On the lattice it is very natural to take the same field 
multiplet, now evaluated at discrete lattice points 



In the corresponding lattice transformations (as defined 
byinD) 

<5C = {eMZ^ + 6-M^)C (82) 

the matrices M and M will be of the same form as in the 
continuum, 



nm 






1 
-V 
-V -1 


0-10' 
-V 

V -1 



(83) 



(84) 



V is subject to the discussion in previous sections. We 
will come back to its locality properties in Section IVEl 



C. Ansatz for the quadratic lattice action 

For the rest of this section we restrict ourselves to a 
quadratic theory. Even in this simple case one obtains a 
nontrivial action with nontrivial lattice derivative oper- 
ators. These operators solve the relation P^ for SUSY 
just as the overlap operator solves it for the chiral sym- 
metry. 

We choose the following matrix K for the quadratic 
action (cf. sect. IIIPp : 





/ -D 


-mb 




-mfc 


-10 


a 





V - m 


\ 


V + m/ 


which means 




S[4>] 1 

a 2 


xUx + V'(V + m/)0 - -FIF 



\ 



(85) 



0n = {Xn,Fn,1pn,'^n) 



n) 



FrubX- (86) 



(again, a sum over the lattice indices is understood). 
Though the symbols of the matrices K suggest them to 
be similar to the objects in a quadratic continuum action 
with W = TOX, they are so far undetermined. (In partic- 
ular the lattice masses and derivatives can be different for 
fermions and bosons.) D, rn;, and nif must be symmet- 
ric and V antisymmetric circulant matrices to guarantee 
hermiticity and translational invariancc. In the contin- 
uum limit we expect the behaviour D — )■ 9^, V ^ 5 and 
rubj — *■ m.1, while at finite lattice spacing these matrices 
are chosen according to relation P^ . This also means 
that the matrix V can be different from the derivative 
operator V in the naive generators M and M . Generali- 
sations of the mentioned ansatz for the quadratic action 
are possible. One could, for example, introduce an addi- 
tional undetermined matrix in the .F^-term, but this is 
not necessary for the solution of the relation in the next 
two sections. In the last section of this chapter we will 
consider the most general ansatz to derive some state- 
ments for the general locality of the solutions. 
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D. First solution for an ultralocal blocking 

To start with the simplest form of the blocking matrix 
a we take its inverse to be 



solution reads 



a{a 





( 0-2 











Jmn 






ao 









ai 




V 





-fll 






(87) 



where all a^ are symmetric circulant matrices of length 
dimension i (and allowed to be zero, cf. Section riVB[) . 
Note that these considerations also include a much more 
general ansatz. This happens because ([57)1 is the same 
as a general blocking ansatz up to a symmetric part as 
as shown in appendix [Fj 

Circulant matrices commute. Using this property, a 
solution of the symmetry relation (|18p for the quadratic 
SUSY action ([55]) is straightforward. In Section IIIDI 
we have already specified the symmetry relation for 
quadratic theories. From the first line of (|27p one can 
read the following equations 

-D + VV + V{mf - rub) 

= - (aaV + ai)n + (aiV + ao)mb (V + mj) (88) 
V — V + mf — rrib 

= - (a2V + ai)TOf, + (aiV + ao) (V + to/). (89) 

Two additional equations can be identified with the 
transposed of these if one reconsiders the symmetric 
or antisymmetric form of the matrices. The field- 
independent second line of ((?7)) vanishes as M and hence 
Ma~^K^ always connect bosons with fermions. The cor- 
responding relation for the generator M induces the same 
set of equations. 

At first we proceed in the same manner as in the deriva- 
tion of the Ginsparg- Wilson-relation: we use an ultralo- 
cal blocking with o^'s diagonal in lattice sites and derive 
a solution for the lattice action in terms of these matri- 
ces. This solution corresponds to the overlap operator, 
that also is a function of the blocking (appearing in the 
chiral case). 

The second equation ([89]) can easily be solved for V -I- 
ruf in terms of mb and V. One gets V and ruf as the 
antisymmetric resp. symmetric part of 



(1 + aiTUb + ao) - {a2mb + ai)V (V + nib) 



X 



(90) 



where 

X = {1 + ainib + ao)^ - ia2mb + ai)'^V'^ . (91) 

The first equation ([55)1 then gives D and the complete 



V = 



TO/ 



(1 -I- ao - a2TO^)V 



X 

(1 + ainib + ao)mb - (a2TOh -|- ai)V- 



X 



-D- 



-V^ 



1 + ao — a2V^ 



(92) 

.(93) 

(94) 



The last part is presented in terms of — D + ml because 
this operator appears in the bosonic sector after integrat- 
ing out F, i. e. in the on shell action. 

As expected, in the limit a^ — > one has D ^ V^ , V — » 

V and ruf —^ nib. At finite lattice spacing, however, these 
operators are nontrivial because each of them contains 
both the derivative operator V and the mass term mb- 

V and ruf are nonsingular because X is positive since V 
is antihermitian; D is nonsingular, if ao and a2 are not 
largely negative and the denominator in (I94|) vanishes. 

Now we have to check the locality of our resulting lat- 
tice action. As explained in section llV) the solution of the 
additional constraint, (|?T)) . leads to the non-local SLAC 
derivative V(p) ~ p. As the operators of the lattice ac- 
tion ([5^ - ([M)l are given in terms of V, one might expect 
that they inherit this locality problem. Therefore, the lo- 
cality properties of the lattice action need to be examined 
carefully. 

Since a is similar to a mass term (diagonal in lattice 
sites) the form of the denominators in ()92p - (|M)) renders 
the appearing operators very similar to massive propa- 
gator. In the continuum similar expressions lead to an 
exponential decay for large distances. On the lattice, 
however, the corresponding behaviour is spoiled by terms 
decaying only algebraically. This is shown in appendix 
l^ni using methods of complex analysis. 



E. Solution with a local action 

Since we insist on the locality of the lattice action, 
more general blocking kernels a must be considered. 
With these it is possible to enforce locality in the SUSY 
lattice action. In our point of view this is a crucial feature 
of the modified symmetry including the blocking kernel 
compared to the naive symmetry without it. 

Allowing now for an arbitrary momentum dependence 
of a one can solve the equations (|55)) and (|55)) for the 
circulant matrices ao, ai and a2. Consequently these 
matrices are dependent on V, D, mj and rub- Given 
V(pfe) — ipk the solution in terms of the matrices of the 
lattice action is 



aoiPk) 
ai{pk) 

a2{pk) 



D 



iVpfc 



D 



-TOf, 



TO/ 



-D 



D 



jy/pk 
m'j-V'^ 



(95) 
(96) 

(97) 
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Now one could use the simplest ultra-local operators 
(without doublers) in the action and read off the cor- 
responding a~^ for a solution of the relation. Additional 
restrictions for these operators appear since singularities 
in the ai must be excluded. A possible singularity ap- 
pears at pfc = if the mass is zero in the theory. To 
avoid this singularity one can use m\y = rnj and D = VV. 
Note that in this way the fermionic and bosonic sector 
are treated in a similar manner. So possible doublers of 
the fermionic sector also appear in the bosonic one, but 
can be removed with the same mass term. The result 
then simplifies to 



ao(Pfc) 



V(V - iVk) 
V2 



TO, 



a2ipk) 



1 + iV/pk 



(98) 



'/ 



and ai{pk) = 0. Just as expected all Ui vanish if the 
derivative operator V is the SLAC derivative and the 
deformed symmetry is reduced to the naive one. 

The simplest solution for the fermionic operators is 
the standard Wilson fermion. With the corresponding 
bosonic operators one finally arrives at 

V(pfe) = - sm{apk) , D = VV (99) 

a 

rribipk) = mf{pk) = to+ -(l - cos(apfe)) . (100) 

In this realisation all possible doublers are removed by 
the Wilson mass terms in the bosonic and fermionic sec- 
tor. 

Such a form for the quadratic lattice action was chosen 
in [23] and other lattice simulations. It has been shown 
that the choice of the same mass and derivative opera- 
tors in the fermionic and bosonic sector leads to a major 
improvement with respect to the lattice supersymmetry. 

Note, however, that a major requirement is still not 
considered. One should not only insinst on the locality 
of the action. To get a well defined representation of 
the symmetry on the lattice also the deformed symmetry 
transformations must be local. 



F. Local supersymmetry 

In order to get a well defined lattice supersymmetry, 
the deformed symmetry operator Mdof has to be local and 
approach the continuum supersymmetry, as explained in 
Section IIIDI The solutions we have obtained so far con- 
tain the SLAC operator, either in the action, eqn.s (19^ - 
([M]) . or in the inverse blocking kernel, eqn.s ([551) - ([55)l . 
This leads to a non-local behaviour in a~^K, which in 
Mdct = M — Ma~^K can enhance or reduce the non- 
locality of M (whereas in the chiral case a non-local 
a^^K immediately induces a non-local Mdcf). Here we 
investigate the conditions under which the locality of the 
action and the deformed symmetry generator A/dot can 
be achieved. We first derive a special solution which ful- 
fils the locality condition (I24p for the action and Mdcf- 



At the end of the section we will argue that this condition 
(and not (P5|) ) is the best one can achieve in the present 
setup. 

For the special solution we consider now a slightly gen- 
eralised form of the deformed supersymmetry and the 
lattice action. Since the corresponding inverse blocking 
kernel has no direct physical implication it is adjusted 
accordingly. 

In the ansatz for the quadratic lattice SUSY action, 
eq.. (|85|) . a general symmetric circulant matrix —/is used 
in the i^^-tcrm instead of the diagonal matrix —1. As a 
deformed supersymmetry generator we take 



Mdcf = 



/ 



V 









lb 








-Vh 


Vf 


-If 















(101) 



If a~^ vanishes in the continuum limit, it follows that 
Mdcf -^ M, and hence hj — > 1 and Vhj -^ dt- 

Demanding that Mdef is a symmetry of the lattice ac- 
tion K, eqn. (|5^. one arrives at the following sets of 
equations: 



(102) 
(103) 
(104) 



h 


TOb 


Vb 


If 


TO/ 


V/ 


V/ 


D 




h 


V 




Vb 


V 




h 


/ ■ 





This can be solved easily, e.g. for the choice ([M|) . (jlOOp 
of straightforward operators in K one obtains / = 1 . In 
that case a natural solution is 



/b=// = l, Vb = V/=V. 



(105) 



With V as defined in (|M|) ). for example, the deformed 
symmetry is ultra-local and obviously approaches the 
continuum supersymmetry. So it seems that there exists 
a local deformed symmetry of the considered local lattice 
action. The reason for the absence of locality problems 
is that the SLAC operator from the blocked lattice sym- 
metry M has not appeared yet. 

The deformed symmetry must, however, not only be a 
local symmetry of the lattice action. A further condition, 
eq. ((3T|) . implies a relation between M and Afdcf- This 
additional condition is a direct consequence of the sym- 
metry relation psp . Since a~^ is so far undetermined it 
is apparently not difficult to satisfy this condition. But 
restrictions of the lattice action, due to e.g. hermiticity, 
impose further constraints on the blocking kernel. In 
the case of supersymmetry these constraint are of great 
importance since a~^ can only connect fermions with 
fermions whereas Mdcf connects fermions and bosons 
with each other. So one has to investigate whether or 
not there exists an a~^ to ensure that this deformed sym- 
metry Mdcf is indeed the result of a blocking procedure 
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and hence fulfils Mdcf = Af (11 - a'^K). As M is not 
invcrtible we use a general ansatz for a^^, namely 



/b2 b[ 

1 I 1 



a{a 








b[ bo 

bi + b'{ 

V -bi+ b'{ 



(106) 



and compare the two sides of Ma ^ 
yielding 



(Af - Mdcf)if- 



i-l){b^~b'() 




/ 

Vb(&i-6'/ 

-V&2-&i -Vb[-bo 

V 

/ (l-/fa)/(V-m/) 0\ 
(Vfc - V)/(V - m/) 


\ 0/ 



(107) 



The lower left block of this matrix equation can be satis- 
fied by an appropriate choice of bo, b\ and 62- The upper 
right block fixes 61 — fe" and in addition implies 



Together with eq. (|104p this yields 



(108) 



It is instructive to recall that these results rely on the 
specific structure of the transformations in the case of 
supersymmetry. They transform fermions into bosons 
and vice versa. By contrast the blocking matrix and the 
action relates fermions with fermions and bosons with 
bosons. In case of other symmetries the transformations, 
a~^, and K would have the same block diagonal struc- 
ture, which makes it easier to find a local solution for K 
and Mdcf. 

One may still wonder whether the current approach 
can be generalised to get a result that satisfies the more 
severe condition (|^5|) for locality. To investigate this 
problem we look at the relation Afdof = —Ma^^K (see 
(PIjl ). where only the symmetric part of a^^ of the block- 
ing appears. To get an exponential decay for both Afdef 
and the action K, also Ma^^ has to fulfil this condition. 



The SLAC-derivative in the generator M, (|85|) . violates 
the locality of this matrix. Taking the most general a^ , 
lESlI, the only entries of Ma^^ are di + d'{ {di - d'{ for 
the second supersymmetry, Mag ) and the same matrix 
multiplied by the SLAC-derivative. So both rfi + rf" and 
V(di -I- d'O must be local. In analogy to / and / V we 
conclude that it is impossible to get an exponential decay 
on the lattice for both matrices. So the condition ([24]) 
and not (1^51) can be fulfilled for both A//dcf and K. 



/V = V 



(109) 



We remind the reader that V is the SLAC derivative ap- 
pearing in the lattice symmetry M due to the additional 
constraint and / and V are the matrices for the F^ and 
the kinetic fermion term in the lattice action, which both 
should be local. 

The SLAC operator is non-local because of its discon- 
tinuity at the boundary of the Brillouin zone. To ensure 
a local V = /V, i.e. to make it periodic in momentum 
space, / and all its derivatives must vanish at the bound- 
aries of the Brillouin zone. Then the behaviour is not 
analytic in momentum space, but a stronger than poly- 
nomial decay is guaranteed (see appendix Fill . With such 
a result for the matrix / all symmetry conditions for the 
action can be fulfilled with 

m/ = TOfc , If = lb = I (110) 

V/ = Vh==v = /v, n = /v^ (111) 

All of these operators of the action satisfy locality 
stronger than polynomial. Note that this solution 
amounts to Afdof = IM. 

The only remaining problem arises since I vanishes 
at the boundaries of the Brillouin zone. This means 
that the on-shell bosonic mass term inherits a divergence 
TOh(p)^/2f (p) at the edge of the Brillouin zone. Hence an 
on-shell problem is expected although the off-shell action 
is local in the sense of condition 



VI. TOWARDS INTERACTING 
SUPERSYMMETRY ON THE LATTICE 



The final goal of our investigations of lattice SUSY is 
to write down a lattice action for an interacting theory. 
According to the usual argument, that is not blocking in- 
spired, a supersymmetric lattice action for the quadratic 
case can be found. The lattice supersymmetry transfor- 
mations are in that case defined as the continuum trans- 
formations with the derivative operator replaced by a 
local lattice derivative. 

For the case of a quadratic action we have given a so- 
lution above, that preserves the modified symmetry and 
is hence guided by the blocking of the symmetry. In this 
case it is also possible to derive a direct solution of the 
blocking transformations, cf. eq. (|T7|) . Therefore it is 
desirable to extend the above results to interacting the- 
ories. 

For actions beyond second order we go back to the orig- 
inal equation (jlSp which provides a systematic relation 
to be fulfilled in order to keep the considered symmetry. 
The solutions are, however, much more complicated than 
in the quadratic case. This can already be observed by 
an analysis of a specific solvable example, namely con- 
stant fields in SUSYQM, which gives nontrivial results 
and displays a new issue: the polynomial nature of the 
action. 
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A. Solution for constant fields in SUSYQM 



In this subsection we work with the same parametri- 
sations for M, M and a^^, see ([55]) and ([57)l . but study 
only constant fields Xn = X s-^d so on. This amounts to 
the zero mode sector of the lattice derivative V, i.e. we 
can replace V — *■ 0. Note that in this approximation F 
is invariant under the naive lattice transformations. Of 
course, locality will not be an issue for this toy model. 

For the action we use the following ansatz: 



- = N[i,i:g{x)~hix,F)] 
a 



(112) 



where N, the number of lattice sites, is a remnant of 
the summation. In view of the continuum limit we have 
restricted ourselves to an F-independent function g cou- 
pling the boson and the fermion, like dW/dx in the con- 
tinuum action (j76p . Likewise, we expect the Fg = dh/dx 
to hold in the continuum limit. 

For such an action the relation p^ becomes a partial 
differential equation in g and h: 



^ dh dh dh ax dg 

Fg~- = -Na,g--Na,g—--- 



(113) 



This indeed approaches Fg — dh/dx in the limit where 
the Oj vanish. 

For finite oq and ax this equation can be solved for 
h for different choices of g. Among these solutions we 
consider those, that consist of a term F^/2 plus terms 
linear in F with an arbitrary ^-dependence, such that 
the auxiliary field F can be integrated out easily. 

The simplest case is g{x) = 0; which should include the 
'kinetic term', that in the zero mode sector degenerates 
to —F"^ /2. Indeed, any function h{F) is a solution of 
eq. (|113p in this case. 

The case g{x) = "i/ resembles an additional mass 
term. The corresponding solution 



Hx.f)^\f^ 



l + aaN 



■mtFx 



1 — aiNiTif 

a„ {1 + aoN)Nm} 
^y (l-ai7VTO/)2 



X 



(114) 



becomes F'^ /2 + irifFx in the limit a^ -^ as required. 
For finite Oj we obtain a bosonic mass different from ruf 
and an additional term proportional to x'^- This action 
could also be gotten from Section IV C[ eqn.s ([55)1 and 
(P^ - ((M)) . since putting V — > there also yields different 
bosonic and fermionic mass and a x^-term surviving from 
D. 

The most interesting case is an interacting theory with 
a truly x-dependent term g{x)- For the lowest possible 
power 



gix) = Ax: 



(115) 



one obtains the common Yukawa interaction. The gen- 
eral solution of (|113p is again restricted by the require- 
ment that for vanishing constants a^ the term h should 



resemble the continuum result F'^/2 + XFy^ /2. One ob- 
tains the non-polynomial solution 

U^ F) = If^ - i±^vF + -o(^ + ^oN) ^2 



x]logil~aiXNx) 



1 l + apA^ ao(l + aoi V) 

aN ajXN^ afXN^ 

which upon expanding in ai becomes polynomial 

Hx,F)^^F^ + ^il + aoN)Fx^ 

+ ^X'Nil + aoN)x^ + Oia,). (117) 

Now it is obvious that the correct continuum limit is ap- 
proached and that this interaction term does not diverge 
in the limit ai — > as one might have expected from 

dni. 

Note that in the zero mode sector there is no failure 
of the Leibniz rule. Hence the direct translations of the 
continuum action are naively supersymmetric fulfilling 
p9)) which amounts to setting ai_2 = (in other words 
using a SUSY preserving blocking as)- The presented 
solutions are deformations thereof. Already in the sim- 
ple case of constant fields these deformations lead to a 
non-polynomial solutions if an interacting theory is con- 
sidered. 

One might expect that the non-polynomial form of the 
action is a generic consequence of the symmetry relation 
psp . not only for supersymmetry. Indeed, similar results 
have been obtained in the context of chiral symmetry 
[rn |. We will discuss this in the next sections as well as 
the circumstances under which a polynomial action like 
(|117p is possible. 



B. Polynomial solutions of the symmetry relation 

As shown in the previous subsection it seems difficult 
to find a polynomial solution of the lattice symmetry re- 
lation (fT8|) beyond second order. In the following we 
argue that a non-polynomial action is indeed the generic 
solution of this relation for an arbitrary linear symmetry. 
Only if special conditions are fulfilled, the series in the 
fields can be truncated. To show this general behaviour 
we generalise the considerations of section IIIDI to inter- 
acting systems. 

We consider a lattice action consisting of polynomials 
up to order R in the fields represented as 



r=l 

where s^''-' contains the rth order in the fields. The purely 
quadratic case R = 2 (and s'^' = 0) has been discussed 



in section |Vl The coefficients K are so far not further 
specified, they can imply a simple multiplication of fields 
at the same lattice point, but are also allowed to contain 
lattice derivatives or to smear the powers of the fields 
over several lattice sites, as long as they obey the correct 



J 
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continuum limit. 

The relation (jlSp is in general a complicated nonlinear 
differential equation coupling derivatives with respect to 
the fields at different lattice points. An expansion in the 
order of the fields using the ansatz (|118p yields 



0(0 



0(/) : 



r=l...R-2 



oi<p' 



r=R-l,B.\ 



O(0 



r=R+1...2R-2\ 



Ma 



.i/,5sW(5s(i) ^2^(2) 



M0- 



M0- 



Ss^"-^ 



{STrM ~ STrAf ) 



Ma- y: ( 

s+t=r+2 



5si-) 5s(*) 52g(r+2) 



Ma- J: 



,5s(^' (5s(*) 



s+t=r+2 



= Ma-^ J: 



Ss^"^ (5s(*) 



(119) 
(120) 

(121) 

(122) 



s+t=r+2 
I 



where we used a short-hand notation without indices. 
These coupled equations can be read as restrictions for 
the ^^^"'^ parametrising s^^' imposed by the symmetry. 
In the case of R=2 only the conditions (|119p . (|120p are 
relevant giving ([?7|) . 

For interacting theories, i? > 2, a set of equations with 
vanishing left hand sides, eqn.s p22p . appears. The very 
difficulty to obtain a polynomial solution is to fulfil these 
equations with a finite number of interacting terms s^"^' [(/)] 
(which in addition shall give the desired theories in the 
continuum). If this turns out to be impossible for some 
order R, this order has to be increased and finally one 
might be forced to non-polynomial interactions. 

As an example we consider the relation of the highest 
order 0{(jP^-^'^), which reads 



<5s(«)<5s(«)- 






(123) 



and can be rewritten in matrix-vector notation as 



= v'^{Ma~^)v with v^ 



9s(fl) 



(124) 



This relation is a severe constraint, because it implies 
that 



Ma' 



{Ma-'f = 



(125) 



within the subspace of lattice fields spanned by the w^. 
If the w^ span the whole space of (f>l^ , the relation is im- 
mediately reduced to the naive symmetry. 

If the w^ do not span the whole space of the fields they 
must be linearly dependent. Then some linear combina- 
tions of the wjj vanish and from the definition (|124[) it is 
clear that the highest part of the action s'^' does not de- 
pend on some particular combinations of fields. On this 
subspace there is no constraint like (|125p . 



So after all it is only possible to get a truncation of 
the action, if (|125p is fulfilled on that subspace of 0's, 
on which the highest term of the action, s'^-^^ depends. 
Keeping translational invariance, it is impossible to have 
g(R) independent of fields at particular lattice points n, 
but s^ ' may be independent of a whole field component 
(0*) of the multiplet. Such a case appears in the case 
of constant fields (see eq. (|117|) ') when oi is set to zero. 
Then the highest term of the action, x^, depends only on 
X, and Ma^^ + {Ma~'^)'^ has no matrix entries for this 
field component. In this way a polynomial solution can 
be achieved. 

We stress that so far we only discussed one necessary 
condition for a truncation of the action down to a poly- 
nomial. The remaining equations (|119p - (|122p need to be 
solved as well. 

The continuum limit of the resulting actions needs 
in general a careful investigation since additional inter- 
action terms must be introduced to solve the relation. 
These terms have no corresponding continuum counter- 
parts and should vanish in the continuum limit. It is easy 
to ensure this in a naive way, where one just introduces 
the appropriate power of the lattice spacing, a, in front of 
the terms to let them vanish when a goes to zero. Notice, 
however, that the additional interaction terms introduce 
new vertices in a perturbative expansion. Because of di- 
vergences (for a — > 0) the corresponding diagrams are 
not necessarily vanishing even though the vertices them- 
selves vanish. So it seems non-trivial to find the correct 
form of the action in the continuum limit. 

The situation gets even more difficult if the truncation 
constraint is not fulfilled. Then one has to accept a non- 
polynomial actions. Beside the feasibility of such actions 
for numerical simulations, one has to analyse fundamen- 
tal aspects of them like renormalisability. Nevertheless, 
these solutions obey the deformed lattice symmetry that 
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approaches the continuum one in the continuum Hmit. 
This hmits the possible actions, e.g. by relating the dif- 
ferent couplings, and might even determine the entire 
form of the action. 

How the discussion specialises to supersymmetry will 
be discussed in the next subsection. 



C. Interacting SUSY 



vanish in the continuum limit, such that the continuum 
action contains only terms of lower order in the fields. 

Again, in this approach we so far have discussed only 
the highest relation. To solve the remaining relations is 
a non-trivial task to be done in further investigations. 

With the help of such additional terms it might be 
possible to fulfil the symmetry relation which then guar- 
antees the (super)symmetry in the continuum limit. 



To illustrate the above rather formal general state- 
ments, we specify them to SUSY theories. For the sake 
of the argument we consider now a higher order action 
of SUSYQM, other SUSY theories will lead to a similar 
situation. 

We mimic the interaction term of the continuum the- 
ory, f/J^X + Px^ (see eq. ([75| with W — x^/2), by the 
following lattice action: 



3 — -'ini,n2,"3 rni'rn2Xn3 



" ^rLi,n2,n3^ niXn2Xn3 



(126) 



with parameters K and K' as in pisp . 

The (super)vector v following from this action reads 



/ ^ni,n2,nWn-iWn-2 ' '^^n-i_,n-2,n^riiXn2 \ 
-^ni,n,n2'YniX'n2 



V 



K„ 



ii'n.Xr, 



(127) 



/ 



It is not hard to see that v takes all values in the con- 
sidered space of lattice fields upon varying the fields 
(Xn, Fn,'ipn, i^n) at all lattice sites. The rationale for this 
is, as explained above, the 'completeness' of the highest 
term in the action s^^\ eq. p26p . if it contains all fields of 
the multiplet. As a consequence, (|125p has to be fulfilled 
in general, but this reduces the symmetry relation (jlSp 
to the naive symmetry. One can easily convince oneself 
that the same argument applies to all continuum-inspired 
actions of the form {R — l)ipipx'^'~'^ + Fx^^^- 

There are two options to proceed at this point. The 
first one is to accept non-polynomial actions like (|116p . 
which induces many new complications. 

The second and rather intricate possibility is to fulfil 
(|125p only on a subspace of the fields. For SUSY this can 
be achieved, e. g., by building the highest order term (in 
this case s^^') out of purely bosonic or purely fermionic 
fields. As Ma^-^ has only entries mixing bosons and 
fermions (|124p is fulfilled easily in this way. We repeat 
the constant field result, where setting ai = and ag ^ 
in eq. pi7p the highest term is of the form x^- In the 
present case with 'smearing' in K^ such a term could be 

repreSenteQ as -'^nl,n2,n3,n4X?^lX?^2X?^3X^4■ 

Of course, such incomplete parts of the actions on their 
own are not giving the desired continuum limit since in 
the continuum SUSY actions all fields are present at any 
given order (i.e. in every s^'"'). The only meaning such 
terms could have is to be lattice artifacts needed to solve 
the symmetry relation (|18p at finite lattice spacing and to 



VII. SUMMARY 

In this work we have systematically approached lattice 
supersymmetry, and have extended the Ginsparg- Wilson 
relation to general linear symmetries and interacting the- 
ories. 

For a general blocking procedure we have derived a 
lattice relation, eq. (HH]), from the continuum symmetry 
of the theory at hand. This relation can be viewed as 
the 'remnant' of the continuum symmetry on the lattice, 
and has to be satisfied by the lattice action. The relation 
also includes potential anomalies and reduces to the well- 
known Ginsparg- Wilson relation for chiral theories. 

As this relation is derived for all possible blockings 
it does in general not comprise a lattice symmetry, but 
rather describes the breaking of the continuum symme- 
try by the blocking procedure. It can be interpreted as 
a deformed symmetry only for those blockings for which 
the symmetry operator in (PT|) is local, as defined in (|24p , 
and tends towards the continuum symmetry operator in 
the continuum limit. These important properties are dis- 
cussed in section III CI These requirements exclude for 
instance the Wilson fermion action as a solution for a 
chiral lattice theory. 

Interestingly enough, the averaging function / defines 
the lattice formulation, but does not appear in the re- 
lation of the lattice action. Instead, / is involved in 
an additional constraint encountered in the derivation of 
the relation. This constraint is of particular importance 
when the symmetry transformations include derivatives, 
which is one of the characteristic features of supersymme- 
try (as opposed to e.g. the algebraic 75-transformations 
in the chiral case). In order to keep the resolution of the 
averaging function as demanded by the lattice cutoff, the 
derivative operator in the lattice transformations must 
be the SLAG derivative or approach the SLAG deriva- 
tive in the continuum limit. Such a non-local operator 
is of course problematic in view of the locality of lattice 
action and symmetry. 

The other ingredient of the blocking transformation, 
the blocking kernel a, appears explicitly in the symmetry 
relation of the lattice action. An a that respects the sym- 
metry leads to a vanishing rhs of this relation and thus 
to a naive symmetry. For chiral theories this is forbidden 
because of vector symmetry. In case of supersymmetry a 
symmetric a can in principle be chosen, but the result- 
ing deformed symmetry would then contain the SLAG 
derivative and hence be non-local. The corresponding in- 
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variant lattice actions inherit this problem. Thus we are 
lead to non-symmetric a's and a non- vanishing rhs of our 
relation to have a chance to obtain local lattice actions 
and local symmetry transformations for supersymmetry. 

In the concrete example of quadratic SUSYQM one 
can find how a improves the locality of the action: An 
ultra-local a^^ introduces propagator- like denominators 
in the lattice operators, which however are not enough to 
compensate for the non-locality of the emerging SLAC 
operator. Considering a~^ as a momentum-dependent 
quantity allows (ultra-) local lattice actions as a solu- 
tion of the symmetry relation. Beside the lattice action 
also the deformed lattice symmetry must be local. To 
get a connection to the blocking procedure, only a re- 
laxed version of locality can be guaranteed (cf. section 
IV Fp . Hence this rather simple example already reveals 
a lot about the interplay of the blocking ingredients, the 
locality of actions and symmetries on the lattice. 

One of the main differences between SUSY and chiral 
theories is that the symmetry of the former acts on non- 
quadratic terms. Correspondingly, the relation for lattice 
SUSY extends beyond second order and couples different 
powers of fields. We have discussed that this generically 
results in non-polynomial actions. This finding is not 
completely surprising as our construction is half way to- 
wards the full quantum effective action: we integrate out 
the original fields and represent the remaining path in- 
tegral in terms of the blocked lattice fields. The non- 
polynomial actions are derived from the continuum the- 
ory in a prescribed way. One might for instance speculate 
that certain cancellations between bosons and fermions 
are still present and help to renormalise the system. We 
have given a necessary condition for a truncation in the 
action. Indeed, in the very special example of constant 
fields in SUSYQM this criterion could be satisfied and 
the action is local. 

We conclude with a short outlook. To date, the 
main problem for performing SUSY lattice simulations 
is the implementation of the symmetry. We hope to 
use the blocking formalism to obtain local interacting 
theories, which are invariant under a local deformed 
symmetry. The corresponding solutions to the symmetry 
relation should either be obtained through reasonable 
approximations or by a resummation. This is work in 
progress. 
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Appendix 

A. Solution for a quadratic action 

The starting point to get the fixed point operator in 
its usual form is the result of the Gaussian integration, 

K = a-af{faf + K)-^fa 

= a[l^f{faf + K)-'fa]. (A.l) 

A direct inversion of f^ af + K is difficult since f^af is, 
unlike K, non-diagonal in frequency space. This property 
comes from the /, that maps a larger space onto a smaller 
one. So f^af acts on a larger space than the diagonal a 
and is in general non-diagonal. Therefore, one first gets 
a closed expression of the inverse of K. Assuming K and 
a to be invertible an expansion in terms of a Neumann 
series yields 

K-' = 5:[/(/^a/ + i^)-V^«]"a-i 

n=0 



x{faf + k)-^f 

oo 

x{faf + K)-'f. (A.2) 
After a resummation the result becomes 



K-' = a-' 



= a-^ + fk-'f. (A.3) 

An inversion of this expression is the know result men- 
tioned in equation (|48p. 

An other way to derive it is done via the introduction 
of auxiliary fields c^ , 



dip e 



'j-fv)a{<p-fv)-^vK'P 



-S[0] _ 



After a Gauf5ian integration (first of the ip then of the a) 
one gets again the desired expression for K 



M' dip /dcre-5^'^^'^+^'"('^-/^)-5'^" ''" 



g-S[0] ^ ^„g-i0(/i?-i/^-Ha-i)- 



(A.4) 



With the same conventions as in section[B]one can easily 
derive the Fourier representation of K. 
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B. Fourier analysis of the additional constraint 

For simplicity we work in a one-dimensional finite vol- 
ume of size L. A general continuum field has the Fourier 
series representation 



with 



V{^) = Yl "^(P?) 



^iPqX 



ViPq) 



L 



dxip{x)e 'P""" , 



(B.5) 



(B. 



with dimensionless wave numbers g (g Z) and pq = -j^. 
The same representation can be applied to the averaging 
function f{an — x). 

We discretise L — Na with an odd number N of lattice 
points and lattice spacing a. Functions on the lattice 
can be parametrised by N independent waves which we 
choose to be in the first Brillouin zone. 



(j){an) 



(N-l)/2 

E 

fc=-(Ar-l)/2 



</.(pfe)e^P'=" 



n=0 



(B.7) 



(B., 



From this relation it is clear that 4>{pk) is periodic in 
Pk, 4>{.Pk) = (t>{Pk + Z27r/a) = (/^{pk+iN) V/ e Z. The same 
transformation is used for the (j)f{an). (For (t)f{pk) k runs 
from -{N - l)/2 to {N - l)/2.) 

In Fourier space the convolution in the averaged field 
of eq. ([2]) becomes a product, 



(f>fiPk) 



oo (N-l) 
n= — oo m=0 



e^iP^-P-)--f{pq)^{pq) 



.27r. 
a 



L Y. f{Pk + l — MPk + l — ). (B.9) 

/ — — OO 



.27r. 
a 



This shows how the averaging projects the Fourier com- 
ponents of if onto the first Brillouin zone. In addition 
one easily observes that the Fourier components of <j)f 
and the lattice fields are determined by /, which means 
that / introduces a cutoff for the lattice momentum if 
f{pq) vanishes for all pq greater than the cutoff. 

Because of the constraint (|56l) . the sum on the rhs 
actually has at most one term. 

The additional constraint ([S^ reads after partial inte- 
gration 

y^'^nmf{O'iTi — x) + dxf{an~x) = Vn,.x. (B.IO) 

m 

Because of the circulant form of V (cf. eq. ([55)1 ) we define 
its Fourier transform as 



{N-l)/2 



V„„. = l Y. V(pfc)e^P'='^ 



(n—m) 



(B.ll) 



(TV-l) 



v(pfc)4i = ^ Y 



^ipkan ipiam 



(B.12) 



m— 0, n— 



where k and / are integer numbers between —{N — l)/2 
and {N — l)/2. For greater values of k the result of the 
Fourier transformation again fulfils V{pk+iN) ~ V(pfc) 
and Ski has the same periodicity (which means it is one 
only for k — I mod N and otherwise zero). 

With these Fourier transforms the constraint becomes 

oo 

E /b<;)[v(P9)-*P,]e'^'^'^"-"^=0, (B.13) 

g— — oo 

which for every individual component pq gives the con- 
straint (l54l). 



C. Doubler problem 

In this appendix we discuss solutions V of the addi- 
tional constraint which are ultra-local. These include the 
(slightly modified) symmetric difference, ((57|) . ([59]) and 
extensions thereof spreading over neighbours at larger 
distances, e.g. (|5T|) . They come with blocking functions 
/ extended over the whole volume, (j60p . or subvolumes, 
depending on how many Fourier components / can have. 
In Section IIV Al we have argued that only the / with 
the highest resolution 0{a) leads to a reasonable, namely 
non-redundant lattice theory. This amounts to V being 
the SLAC operator. 

Although the use of these ultra-local V's is question- 
able, we nevertheless want to point out that doublers in 
them can be removed. For supersymmetric theories it is 
natural that also bosons have a doubling problem [23] , so 
we investigate the kinetic operator in the bosonic sector, 
too. 

Consider the solutions ((92)) - ((94)) . The operator V in 
them has doublers at the edge of the Brillouin zone, un- 
less V is the SLAC-operator, and V inherits them via 
the solution ((5^ . 

To resolve the doubler problem we make use of the fact, 
that the matrix irib in our ansatz ((55)) and in ()92p - ((M)) 
is only restricted to be symmetric and circulant. In the 
spirit of the Wilson term, we will now replace mb by 

V 
{'>TH,)nm — mdnm + 7;— (2^,i„i — 5n+l.m — ^ri-l,m) (C14) 

za 

with r a dimensionless parameter. Such a momentum- 
dependent correction will then also occur in m/ through 
our solution^ and all - bosonic and fermionic - doublers 
are removed as long as the parameter 02 remains small. 



fc=_(Af-l)/2 



'^ In our setting, where m/(mi,), it is the simplest to change mt,. 
One could also have solved mi,{mf) and then change ruf. 
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The dispersion relations of mi, and V are 

T 

mb{pk) = m + - {1 - cos{apk)) 
a 

. {N-l)/2 



V{pk) = - Vl Q sa\{lapk] 



(C.15) 
(C.16) 



1=1 



We consider massless fields to = 0. Then because of 
TOb(O) = and V(0) = it is clear that the bilinears 



-x{-U + ml)x, V5(V 



mf)'ip 



(C.17) 



according to (P^ - ([M)) have zeroes at vanishing p as they 
should. At the edge of the Brillouin zone, p = vr/a, they 
would vanish as well because of V(7r/a) = 0, if there 
were no r-corrections. The finite value mb{T:/a) = 2r/a, 
however, removes the doublers since 



i-D + ml)ip=-) ^ - 
a \ a 



1 



1 + ao 
1 



(V+TO/)(p=- = — -— ^-— 

a a l + zrai/a + flo 



(C.18) 
(C.19) 



These values are non-zero and scale depending on the 
behaviour of the a.;. 

What remains to be shown is that there are no other 
zeroes at intermediate values of p. For that we specialise 
to the solution of V ranging to ncxt-to-nearest neigh- 
bours with ci^2 from eqn.s (|62p . ((63)) and plot in Fig. [3] 
the operators — D + rn^ and V-Hto,^ in momentum space. 
As — D-f TO^ and m/ are real symmetric (also after includ- 
ing r), their Fourier transforms are real, whereas that of 
the real antisymmetric V is purely imaginary. In order 
to seek for zeroes, we therefore plot (— D -I- m^)(p) and 



|(V -I- ■mf){p)\ — J^'^(p) + m'i(p), both as functions of 
p. As one can see from that figure, the doublers are re- 
moved and no additional zeroes appear, unless 02 is too 
big. 

D. Alternative derivation of the symmetry relation 

Let us discuss how the symmetry relation emerges in 
the setting of section IIVBI with currents. Here we take 
the point of view that the Wilsonian action 5* is defined 
to give the same generating functional Zf[J] in (j67p . 



Zf[J] 



1 



-si<p]+j<p 



for all currents J as the classical action in dSS]! , . 



(D.20) 



ZfiJ] = Jf J ^'P^ 



-SdM+j0/M 



This yields the same expectation values for all observ- 
ables composed out of the blocked fields (pf . 



Now performing a symmetry variation ip —> if + eMip 
on the classical fields ip and changing the integration vari- 
able in (jD.21|) accordingly, one gets to 0(e): 

The second term can be rewritten as eJM(J3f[^] using the 
additional constraint. Then it is just the change of Zf[J] 
under J ^ J + eJM, 

eSTiMZf [J] ^ Zf[J + eJM] - Zf [J] (D.23) 



This can be computed in the representation (jD.20p . 
where again we change the integration variable — > 
(j)+eM(j). This yields a STrM part, the variation of S and 
a bilinear term in J from the variation of the prefactor: 



= e<^ STrAf - STrM + SS - JMa 



-'j}Zf[J] (D.24) 



The currents in the last term can be rewritten as 0- 
derivatives of exp(J0) in Zf[J] and, by partial integra- 
tion, of exp(— 5" [(/)]): 



= e 



/ d(/ie"^f'^l+'^'^|sTrM - STrM + SS 

]Zf[J] (D.25) 



,SW^Ma~'^e'S 



Demanding this equality for all currents J, the curly 
bracket in the integrand has to vanish giving (fT7|) and 
thus our symmetry relation (jlSp . 

Notice that we have used the same ingredients as in 
the original derivation of the relation, including the ad- 
ditional constraint and the inverse blocking kernel a~^, 
but we have not relied on the definition (|68p of the Wilso- 
nian action based on a itself. 



E. Singular a 



(D.21) with 



The results of section HV B I extend to a's with some di- 
verging matrix elements amn- The corresponding eigen- 
modes a(/)sing = Aging 0sing with Aging -^ 00 are fixed: 
0sing — (t'sing,{[v]- The remaining modes in (j) still un- 
dergo the smearing procedure. 

The other extremal case is PsingQ^~^^sing — ^ 00, where 
Psing is the projection operator on the singular part 
of a~^. Still (|M|) is vahd but the definitions of the 
Wilsonian action (p5)l and the normalisation (p5)l de- 
velop singularities. If we want to keep these definition 
we could simply change the generating functional Zf[J] 
to Zf[{l — Psing) J] which generates correlation functions 

{4'f,reg- ■ •4>f,rcg)- Here (j)f, rcg = Prcg4>f with P^g = 

(1 — Psing)- Then the above derivations are unaltered 



— J rcgO^-'reg aUQ O^reg 



P, 



reg<- 



Q^ -^rcg ■ 

With this procedure we have removed the (^sing-modes 
from our theory. In particular the definition of the Wilso- 
nian action (j68p with aicg leads to S'[0] = S[(f>i-cg]. 
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a2(_n+m2)(pfc) 
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FIG. 3: Dispersion relations for the bosonic and fermionic bilinears with ultra-local V and Wilson mass in mi,, cf. equations 
(|C.15|l and (|C.16[) . on an A^ = 101 lattice. In the top row we plot the bosonic — D + mj and in the bottom row the fermionic 
|V + mf\, according to the solutions (|92|l - (|94p . as a function of the momentum p in the first Brillouin zone. On the left: r — 
with doublers visible at the edges. In the middle: the doublers are removed by r = 1 and (ao,ai,a2) — (0.1,0.1,0.3). On the 
right: additional zeroes in the bosonic sector occur for the choice r = 1 and (ao,ai,a2) = (0.1,0.1,2.0). 



F. General blocking matrices ac 

The restrictions of the blocking matrix originate from 
the hermiticity of the lattice action. In addition, the most 
general blocking matrix ac should connect fermionic and 
bosonic fields only among each other. Therefore, the form 
of its inverse is restricted to be 



a [Ur 



fb2 b[ 





b[ bo 








6i + 


V 


-foi + 6'/ 



b'l 



(F.26) 



where 6" has to be antisymmetric, whereas all other ma- 
trices are symmetric. In order to get translation invari- 
ance all matrices must be circulant. 

The most general symmetric matrix, fulfilling (|20p with 
the symmetry operators M and M defined by (|55|) and 



a(as^)m« = 



v-i< 


di 


di 


V< 











di + 



d'l 






-di + d'l 




(F.27) 



where d'{ is an antisymmetric and di a symmetric circu- 
lant matrix. Such a symmetric part, a~^ , can always be 
added to a'^ without changing the relation (ITT)) . So the 



matrix elements of 

a"^=ae^+as^ (F.28) 

can be reduced without any loss of generality concerning 
the modified symmetry using the freedom of choosing d'{ 
and di . We choose them to be 



di 



-hi and d'l = -h'{ 



(F.29) 



and define 



ao:=6o-V5'/, ai:= 61 + 61, a2 := b2-'^~'^b'{ , (F.30) 

where all ai are now symmetric circulant matrices. The 
resulting matrix with reduced matrix elements compared 
to Uq is just the one given in equation (j87p . Since a^ is 
not important for the deformed symmetry of the action 
the result is not changed compared to a^ . However, 
the deformed symmetry iWdef depends on this part. So 
if one wants to investigate the locality of Mdef the most 
general blocking matrix a^ instead of the a~^ must be 
considered. 

From these considerations it is clear that all entries of 
a~^ can in principle be zero after splitting off a^ . Then 
the blocking matrix would consist only of a symmetric 
part as- This was done in |23| where a~^ = a~^ with 
d" = was used. The lattice symmetry is in that case 
reduced to the naive one. 



G. Decay of SUSY lattice operators 

In this appendix we discuss the locality properties of 
the operators (I92p - (|94p with momentum- independent a^, 
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which due to the denominators resemble massive propa- 
gators. Nonetheless, we will demonstrate that the posi- 
tion space representations of these operators do not decay 
exponentially on the lattice. 

Let us, as common in the discussion of locality, con- 
sider the limit of large lattices. That is, we replace V in 
eqn. ^ hy pk for fc = -(A^- 1)/2 . . . (TV- 1)/2 and take 
N to infinity. Then for the position space representation 
of the (circulant) fermionic operator V 



{N-l)/2 



i ^ , ^'^ 



n—m. — l ; 



(G.31) 



l=-{N-l)/2 



cf. ([55)1 , the discrete Fourier transformation turns into a 
Fourier integral 



{N-l)/2 



Cl 



lim 



- E 



V(pfc)e*P'= 



fc=-(Af-l)/2 



TT /a 



ipla 



dp V (p) e 

Tr/a 

where p is a continuous momentum and 
(1 + ao - a2ml)p 



V(p) 



(1 + aiTTib + a^y + {a2mb + ai^p^ 



(G.32) 



(G.33) 



from the solution 

In the limit a — > and with p- independent a^, the 
Fourier transform gives an exponential decay because 



dp 



P 



Vi + ^2^^ 



^ipla 



exp 



-sign(0^).(G.34) 



\V2\ 



At finite lattices the Fourier transform (|G.32[) can be 
computed in the complexified p-space. The Fourier in- 
tegral consists of contributions from poles where the de- 
nominator of V(p) vanishes and an additional part from 
the corresponding complex contours that encloses the 
poles. 

Because of the regularity of V(p) on the real axis, the 
poles appear away from it. For positive I we enclose 
the poles with positive imaginary parts by adding paths 
p — ±7r/a + i(, (^ £ [0,oo). The residue at each pole 
is some number coming from V, which we assume to be 
non-zero, times a term decaying exponentially in / from 
the Fourier factor in (|G.32[) . For the same reason there 
is no contribution from the asymptotic contour at large 
Im p parallel to the real interval [—n/a, ir/a]. 

The two new contours contribute 



± / dCV(p)e'P"' 
/o 



-±7r/a+z(^ 



(G.35) 



which turn into Laplace transforms 

dC [V(7r/a + iC) - V(-7r/a + <)] e""^'" . (G.36) 



For large spatial distances I only the value of V at p 
7r/a remains, namely 



[V(Va)-V(-Va)] 



l/la 



1 



dC{l-Cla.) = V(7r/a) — 



(G.37) 
where we used that V is odd. Hence the additional con- 
tours give algebraic corrections to V in position space, 
unless V vanishes at the boundary of the Brillouin zone 
p = n/a. Of course the problem of algebraic tails is ab- 
sent, if V is ultra-local. Here, however, we have consid- 
ered V as a function of the SLAG operator V according 
to the solution ([5^ . The way out, which is followed in 
the body of the paper, section FV El is to consider general 
momentum-dependent blocking kernels. 

We remark that the discussion is similar for negative 
I (where one encloses the poles with negative imaginary 
parts). Furthermore, the algebraic corrections do not 
appear in m/ and — D -I- m^, eqn.s ([M|) . ([M)) . because 
these are even functions in p. 



H. Locality condition 

All of the following one dimensional considerations can 
be extended easily to higher dimensions. The considered 
lattice operators O are, because of translational invari- 
ance, circulant matrices, 

0,nn = On-rn = F{a{n - to)) . (H.38) 

The slightly modified condition for locality demands that 
F decays faster than any polynomial. That means 



''F(x)| < oo \/r e¥l, x,y G a¥i . 



(H.39) 



If the Fourier transform of i^(2;), f{p), and its derivatives 
have no singularities and fulfil periodic boundary condi- 
tions at edge of the BZ the following estimation can be 
made 



\x^Fix)\ 



BZ 



{d;fip))e"' 



< 



\d;f{p)\<Cr<^. (H.40) 



BZ 



Consider now a non-local operator similar to the SLAG 
derivative. This non-local operator should have no sin- 
gularities within the BZ for all of its derivatives. The 
boundary conditions are, however, not periodic. Accord- 
ing to the discussion of the locality of K and Afdct it 
should support the modified locality after a multiplica- 
tion with a local operator. In view of the above argument 
the boundary conditions must hence be enforced by this 
local operator without spoiling the differentiability of /. 
Its representation in Fourier space, I{p), must therefore 
be a function that vanishes together with all its deriva- 
tives at the edge of the BZ. In addition no singularities 



should appear within the BZ for any of its derivatives. 
One function that fulfils these requirements is 



Hp) 



lo 



\p\ < e -.v. ^ ^ 
with e < — . 

bl > e " 



(H.41) 



It is clear that /(p) can not be analytic since any analytic 
function that vanishes with all its derivatives at a specific 
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point must be identical to zero. So the common definition 
of locality in terms of analyticity in momentum space can 
not be satisfied. 
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